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Abstract 

Gauge invariant adiabatic perturbations in the cosmology of a scalar field (p{x, t) gravitationally 
coupled to a perfect fluid (radiation or matter) are investigated. The potential of the scalar 
field is of exponential form. This model recently has received much attention in the context 
of 'quintessence' due to new observations of high redshift supernovae indicating an accelerating 
universe. The perturbations are treated in a gauge invariant formalism to systematically avoid the 
residual gauge modes in synchronous gauge. Two important attractor solutions of the unperturbed 
equations are investigated: the first corresponds to the case where the energy density of the 
unperturbed field p(^^-) dominates the cosmology after a short time, the second to a solution 
with a time independent ratio /5(<^)/phyd- The main question addressed in this paper is if the 
inhomogenous fluctuations can be ignored as often argued or not. Indeed, in some important 
cases the perturbations do not die out with time meaning that the ratio of energy density of 
the inhomogenous fluctuations to the unperturbed part is constant after a while. Consequences 
for microwave background radiation and structure formation in the matter dominated era are 
considered in a straightforward manner and it is concluded that for similar initial values the 
anisotropy of the background radiation should be stronger in this model than just for radiation. 
Furthermore, structure formation should have happened slower than in standard cosmology. 
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2 THE "BACKGROUND" SOLUTIONS 



1 Introduction 



In this paper gauge invariant perturbations in the context of quintessence are investigated. Since 
Perlmutter et al. and Riess et al. [0,|2| have pubhshed their results that the universe is accelerating 
instead of decelerating, much work has been done to find the reason or at least a correct kinematic 
description for this behaviour. One of the best candidates is the so-called quintessence matter- 
component that effectively can be described by purely gravitational coupling of a scalar field to the 
standard matter (radiation or non-relativistic matter) There are many different suggestions for 
the form of the quintessence-potential (see for instance One of the first ones is the exponential 

potential as it was originally discussed by Wetterich [Q-^ and Ratra and Peebles This potential 
will be investigated in the form: 



(1) 



In this definition M is the reduced planck mass (M^ = Mp/167r) and a the only free model parameter. 
The amplitude V can be adjusted by a shift of the field <p and therefore does not play an important 
role. This potential is not compatible with observation [11| which could be changed in late cosmology, 
but it still serves as a suitable toy model. The naturalness of such models has been discussed 
by Hebecker et al. Adiabatic Perturbations of this model in synchronous gauge were already 
discussed by Ferreira and Joyce |13, 14| and Amendola |T^, 16] but the influence of inhomogenous 
perturbations onto the CMBR has been ignored. 



2 The "background" solutions 

Prom the Einstein equations. Energy-momentum conservation and the Klein-Gordon equation the 
three field equations for just time dependent quantities can be obtained for Robertson- Walker back- 
ground 1^: 

(2) ^(^t) + 3H{tMt) + ^^ = 0, 

(3) ^ (V{^) + ^^\t) + p(t)) = H\t) , 

(4) pit) + ^Hit)(pit)+pit))=0. 



This work is in an inflationary context, so the density parameter = ptot/pc is set to 1; an 
extension to the other cases would be straightforward. An equation of state is needed for the perfect 
fluid component and therefore the convention 

(5) 3 {p + p) = np . 

will be used such that n = 4 gives a radiation dominated and n = 3 a matter dominated universe. 
The nonlinear equations admit two well-known "attractor solutions". 

On the one hand this is the so-called "tracker solution" (which will be referred to as attractor 
solution I): 

2M t , ^ 2 

6 m = ^0 + In - , H{t) = — , 

a to nt 

= ^ (i) , = ^^^^ (*) • 
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The name arises because the time dependence of the scalar field energy density tracks the behaviour 
of the perfect fluid. This solution exists (meaning it is reached after a short period for virtually all 
initial conditions) for a} > n/2 and there is no induced cosmological constant and consequently it 
does not allow for an accelerated expansion of the universe. 

On the other hand there is the solution (which will be called attractor solution II): 

2M t ^ 1 

^{t) = ^o + In-, Hit) = — 

a to a'^t 

( t\~^ .r, . / 3 \ 1 

P{t) = Po - , V{t) 



toy ' ^ ' a2 J ' 

It corresponds to the case where p » V and exists for < 3, but in the regime 3/2 < a < \/3 
for n = 3 the attractor solution I is dominant. The energy density part of a dynamical cosmological 
constant in this case is non- vanishing: 

/ 3 6 

(7) PA ■■ 



Nucleosynthesis gives a constraint for a during that period of a > \/22 ~ 4.69 Q or a > 4.5 



5.5 1 14]. Therefore, a scenario seems to be appealing where a is not fixed but decreasing with 
cosmic time. Then during radiation domination attractor solution I was valid with a ~ 5 at the 
time of nucleosynthesis. At the beginning of matter domination, attractor solution I still was valid, 
but recently the value of a decreased below a < \/3/2 and the universe "switched over" towards 
attractor solution II. This could explain the acceleration of the universe today through the effective 
cosmological constant. 



3 The gauge invariant perturbation equations 



The problem of residual gauge modes is known to cosmologists for a long time (see for instance |jT^ ) , 
but has attracted general attention in the early eighties when Brandenberger et al. showed that 
indeed many of the generally agreed results on perturbation growth were not correct due to residual 
gauge modes. These cannot be eliminated in general within synchronous gauge. A sophisticated 
solution to this problem already had been offered by Bardeen a few years earlier. He first took 
the approach to eliminate the gauge modes by making a transformation from the original variables 
to new ones that explicitly do not change under residual general coordinate transformations. This 



approach is taken in this paper, following Mukhanov et al. |2C]. The most general scalar metric 
perturbations are parameterised by the four variables (j), ip, B and E within the line element 



(8) 



1 + 20) dt^ - R{t)2B\idx'dt - (1 - 2^) % + R^{t)2E\ 



dx^dx^ 



Here the symbol | denotes the covariant derivative on the hypersurface of constant time, R is the 
scale factor and t denotes cosmological time. Gauge invariant linear combinations of the so-called 
metric potentials <j) and ifj can be constructed by 



(9) 

(10) 



B-RE]R 



. d r 

^ = ^ - r(b - RE 
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3 THE GAUGE INVARIANT PERTURBATION EQUATIONS 



Gauge invariant scalar perturbation quantities for the perturbations Xi dPi dp of the field, the fluid 
energy density and the fluid pressure respectively, can be constructed analogously: 



(11) 
(12) 

(13) 



ip{t,x) := ip{t) +xit,x) , 
p{t,x) := p{t) +6e{t,x) , 
p{t,x) := p{t) + 6p{t,x) , 



X 



(gi) 



X 



+ R(p{b - RE) , 



5e(si) =Se + R'p[B- REj , 



6p 



(gi) 



6p + Rp(B- RE 



The gauge invariant perturbation for a three vector like peculiar velocity becomes: 



(14) 



Ui{t, x) := Ui{t) + 6ui{t, x) 



6u. 



(gi) 



5ui + R[B - RE 



By using these quantities, it is possible to construct the gauge invariant Einstein equations: 

i?2 



-miR^ + R-^] + 



(15) 

(16) [R^ + R^ 
(17) 



4M2 



6T 



.0(gi) 



R ™o(gi) 



4M2 



{2RR + $ + + i?^^' + mm -K^ + ^AD 



4M2 ^' 



In this equation fC denotes the spacial curvature in the Robertson- Walker metric. In the same manner 
a perturbed gauge invariant Energy-momentum tensor tensor can be derived. Since they couple only 
gravitationally, it is just the sum of fluid and scalar field components: 



(18) 
(19) 
(20) 



^-j.0(gi) 



6t: 



.5p(gi) + (^2^ _ ^^(gi) + v,^ x^^'M 



For further discussions the curvature /C is set to zero and it is possible to set $ = 'J/ because all 
spacial off-diagonal elements vanish in the perturbed Energy-momentum tensor. From ( p!5|) and (|lj 
one gets an equation that gives a relation between $ and x^^^^ • 



(21) 



n 



4M2 



n — Q - 2^ 6 — n. 

V3 $ + 



—^x^'^'^-'^v^^x 



(gi) 



This equation does not depend on ^e^'^, 5pfe') or (5n-^'\ This fact is quite remarkable, meaning that 
all dynamical information already is encoded in the two variables $ and x^^''' and their derivatives. 
The fluctuation in the hydrodynamical part still can be obtained from the (O,0)-component of the 
Einstein equations. 

In addition to this equation one needs the perturbation of the Klein-Gordon equation 



(22) 







in order to solve for the dynamics |21|: 
(23) 



^(gi) + 3i7^(gi) _ -^Ax^') + V,^^ x^^^ - m + 2V,^^ = 0. 



5 



4 Discussion 

4.1 Attractor solution I 



Equations (|2lD and (^) are getting fourier-transformed (denoted by " ) and attractor-solution I is 
inserted. The equations are rescaled into logarithmic time r := ln(t/to) and the equations are made 
dimensionless by the substitution x ■= x/^ ■ The two equations are of second order in time and it 
is possible to rewrite them as a system of first order matrix ODE's: 



(24) 









1 





o\ 


' X 






fi 


h 


h 


U 


~(gi)' 

X 















1 


!> 


V ^' ) 






k 




h) 





with < 



fi 

h 

f7 



2n-12 



24-4n 
an 

6 — 71 



n—p 
~3a^ 



+ 



{3-n)k^ ^(2-±h 



n— 6 



a 

6— n 
6a 

-n-2 



and A;r := i?o^o • The prime denotes the derivative with respect to logarithmic time. The per- 
turbation of the hydro dynamical part can be obtained from the (O,0)-component of the Einstein 
equations: 



(25) 
6e 



(gi) 



12a^ — 6n 



n^a x(gi) 
X ■ 



2 

n a 



\2a? — 6n 



^(gi)' 
X + 



n 



6a^ — 3n (6a^ — •in)k^ 



2na^ 
2^ 



n 



It is common practice in cosmology to use the ensemble hypothesis for the universe (see for instance 
[22 1 or ||2^) and so the (4 x 4)-matrix needs to get rewritten as a (10 x 10)-matrix for the correlation 
functions: 



(26) 



d_ 







( ° 
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^ {XkiXf:^)^ 
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fl 
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{Xk^^k-^ 






















1 
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fs 
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fl 
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fl 


fr 


f2 + f8 
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/ 







The energy density of the scalar field in second order can be separated into an unperturbed and a 
perturbed part: 



(27) 



hp^ + + \x' + y{^) + y'mx + )^"{^)x^ 



"f^kin.flu 



f'pot 



' pot, flue 
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4 DISCUSSION 



In terms of correlation functions, this takes the form: 



(28) 



/'kfn,fiuc(^' 



^[x{x,t)x{x,t) + ■^dix{x,t)dix{x, t)^ 



+ 00 



(2^ 



+ 00 



{xx){q,t) - -^Qiih - qi){xx){q,t) 



(2^ 



(29) 



and 

/'pot,fluc(^' 



Air 
(2^ 



47r 
(2^ 



[xx){q,t)q^ + ■^{xx){q,t)q'^ 



-V"{^){xx){q,t)q^dq . 



dq 



The assumption of gaussian random variables and therefore uncorrelated fourier modes was used. In 
order to do a numerical integration of the ODE's ([2^ ) one needs the initial conditions for every fourier 
mode. It is a prediction of inflation that the spectrum will be scale invariant (see for instance |23| 
or ||2^), therefore the amplitude of every correlation function A;- mode should be the same at the 
time the modes are crossing the horizon. For attractor solution I this does happen at 



(30) 



TH = In 



k 



kR2TT 1 



The question how to choose the overall normalization of the spectra remains, and indeed one does 
not find a theoretical answer, because the underlying quantum theory is unknown. In this paper the 
main question is, how the fluctuations will qualitatively evolve in time. Therefore, the estimate from 
structure formation that the fluctuations of the radiation part compared to the background at the 
time of decoupling must have been of the order 10"'^ should be sufficient |^2|. For the sum of the 
correlation functions in A;-space this corresponds to: 



(31) 



(5estr(^)<5estr(^) 



PstT 



1 Att 



dq q^ {5estr{q)5estr{q)) ~ 10 



This is taken as a normalization for the spectra. The fluctuations in non-relativistic matter and 
radiation at the time of decoupling are related by the simple equation: 



(32) 



'5emat(^)5emat(x) 



5estr{x)5estr{x) 

PstT 



The result of the numerical integrations is shown in figure |T] for a = 2.5. 100 different logarithmically 
distributed fc-modes have been integrated from tq = to r = 1 with a fourth-order Runge-Kutta- 
method using a stepsize of At = 0.01 and a maximum fc-value that corresponds to the horizon at 
Tq. The scale invariance has been implemented numerically. 

It can be seen in the plots that during the radiation dominated epoch (n=4) the perturbation in 
the kinetic energy density settles down to a constant value compared to the background whereas the 
perturbation of the potential energy dies out exponentially. This is an interesting result because it 
shows that during the radiation dominated epoch (if attractor solution I was valid) the perturbation 
of the scalar field cannot be neglected as often argued. For other values of a the qualitative behaviour 



4.1 Attractor solution I 
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logarithmic scaling 
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Figure 1: Ratios of kinetic and potential scalar field energy densities for attractor solution I with 
a = 2.5 and n = 4(top), n = 3(bottom). 



a 


1.414 


2.0 


2.5 


4.0 


5.0 


10.0 




1.41 • 10-" 


1.76 ■ 10-6 


8.28 • 10-'' 


9.32 • 10-5 


2.54 ■ 10-* 


4.66 ■ 10-3 



Table 1: Ratios of kinetic scalar field energy densities for attractor solution I. 



is the same, but the value that the kinetic energy density ratio settles down to highly depends 
on a, as shown in table |^. For a value of a = 5 as favored by nucleosynthesis one already gets 
('^kfnfiuc)//'kfn ~ '^■^^ ' 10~^- The direct comparison of the fluctuation in the energy density of the 
scalar field perturbation and the radiation perturbation shows that this ratio is nearly a-independent 
and after some time takes on a value of {Sp\^l{x)5p^l{x)) / (^de{x)6e{x)'j 1.7. This means that the 
fluctuation in the kinetic part of the scalar field energy density is indeed greater than in the radiation 
energy density. Additionally, it is remarkable that after some time the scaling behaviour is the same 
for both (like for the background energy densities). The consequences for CMBR will be discussed 
later. 

During the matter dominated epoch (n=3) a different picture arises from the integrations: both 
kinetic and potential energy densities of the scalar field perturbation decay exponentially compared 
to the background. This means that even if the perturbation of the scalar field has been important 
during the radiation dominated epoch, in the matter dominated epoch it became negligible. Still 
there will be an interesting effect on the growing perturbation modes in the matter density. 
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4 DISCUSSION 
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Figure 2: Ratios of kinetic and potential scalar field energy densities for attractor solution II with 
a = 1.2(top) and a = 0.8(bottom). 



4.2 Attractor solution II 

The gauge invariant perturbation equations can be written in an analogous form like (|2^ ) but with 
coefficients: 




Furthermore, the equation for the time of horizon crossing is changing compared to (^): 



After these two alterations the equations can be integrated analogously and the energy density of 
the scalar field perturbations can be calculated. Attractor solution II becomes valid for a < 3/3 and 
there is a qualitative difference in the behaviour for a > 1 and o < 1 (for a = 1 the time dependence 
of the coefficients and therefore the fc-dependence vanishes and this case is rather uninteresting). 
These two different cases are plotted exemplary for a = 1.2 and a = 0.8 in figure ||. One sees that 
the ratios of the energy densities for a > 1 decrease exponentially, while for a < 1 this does happen 
just for the kinetic energy density ratio. The perturbation in the potential energy density of the 
scalar field settles down to a constant value. Naively, this might explain a part of the missing dark 
matter, as being stored in the local fiuctuations of the scalar field. Attractor solution II could explain 
today's accelerated universe, as already mentioned. But there then must have been a transition from 
attractor solution I in recent times, meaning a long time after the beginning of matter domination. 



4.3 Possible Consequences for CMBR and structure formation 
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Figure 3: Comparison of the spectra of standard cosmology and the quintessential model with at- 
tractor solution I (normalized, n = 4, r = 4, a = 2.5) 



During matter domination the ratio of the energy densities in the scalar field dies out exponentially. 
Therefore, the value of the potential energy density perturbation of the scalar field should be very 
small when attractor solution II becomes valid, and consequently the energy density stored in the 
perturbations is negligible. 

4.3 Possible Consequences for CMBR and structure formation 

The last scattering of the CMBR occurred approximately at decoupling. Therefore it should be 
generally possible to detect the effect of the scalar field perturbation from the fluctuation in the 
radiation part. It is instructive to directly compare standard cosmology with perfect fluid plus 
scalar field. The main problem is obvious: the initial normalization of the fluctuations has been 
chosen to approximately fit the fluctuations of the CMBR (^), so it would only be possible to 
detect a difference between the two models in the energy density if the qualitative behaviour of the 
radiation fluctuation ratio would be different. Therefore, also the case of standard cosmology has 



been examined in the "language" of the differential equations ( pi] ) and correlation functions (26). 
The result is that indeed one cannot recognize a qualitative difference between this and the case 
with scalar field and attractor solution I. This is rather discouraging but there is still the hope that 
there might be a difference in the A;-spectra. To investigate this possibility the condition of scale 
invariance is given up at this time and different initial spectra are considered. Naive power spectra 
of the form 

(35) Ps,s,{k,T)^k~J^'e'^ 

are used, the free parameter is r. They are compared to the respective spectra without scalar field 
in figure ^. The difference between the two models is smaller than the difference between different 
initial spectra and is getting even smaller the earlier the modes have crossed the horizon to the 
inside (corresponding to greater values of k). The figure shows this for r = 4, for increasing time the 
differences are also getting smaller. Therefore, one cannot distinguish the spectra of the fluctuation 
in the radiation density of the two models from these naive considerations. 
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5 CONCLUSIONS AND PROSPECTS 



a 


1.414 


2.0 


2.5 


4.0 


5.0 


10.0 


Std.-Cosmology(a oo) 


y 


0.5478 


0.9996 


1.1282 


1.2558 


1.2838 


1.3204 


1.3324 



Table 2: Fitted growth exponent y for matter density in matter dominated universe(n = 3), accuracy 
~ 10-3. 



One also wants to consider the influence that the existence of the scalar field and its fluctuations 
have on structure formation. The root for structure formation is that during matter domination 
(n = 3) the fluctuations of the matter density are growing compared to the background. This is 
already known from standard cosmology, and it also can be seen in the integrations of this paper for 
99 = 0, the growth of the correlation functions (5emat(2^)<^emat(2^))//5^ is proportional to exp(?/T) with 
y = 4/3 (for 6etnat/p and normal time this is the familiar t^/^ behaviour of the growing modes). In 
the model with scalar fleld and attractor solution I (n = 3) the growth exponent y is changed. This 
change is a-dependent as shown in table |2[ Therefore, the growth of the matter density fluctuations 
is smaller than in standard cosmology for all values of a and which should lie in the region between 
3/2 and ~ 5 for structure formation as discussed earlier. 



5 Conclusions and Prospects 

There are two important cases where inhomogenous fluctuations cannot be ignored: flrst for attractor 
solution II and a < 1 (which could play a role in late cosmology) and second in a radiation dominated 
universe plus quintessence with attractor solution I. The energy density ratio of kinetic scalar fleld 
fluctuations compared to the background settles down to a constant value, this numerical value 
highly depends on the model parameter a. As a consequence, the inhomogenous fluctuation in the 
radiation should have been stronger at the time of last scattering than without scalar field. Sadly, 
the initial normalization of the spectra is unknown and therefore, it is not possible to decide between 
the two models from such simple considerations. The form of the spectra also does not help very 
much. It might be useful to derive the spherical multipole expansion of the fluctuations on the one 
hand and to consider secondary anisotropics on the other. But this seems to make sense only if a 
more realistic description for matter and radiation is taken that pays more attention to the different 



particle species, see for instance [|l4| , |25| ,Pq| and a potential which is fully compatible with observation. 



Furthermore, it might be interesting to consider entropic perturbations in addition to the adiabatic 



perturbations. Appropriate Ansatze can be found in [20|. 

The growth exponent of the growing modes for structure formation in the matter dominated 
universe is changed by the coupling to scalar fleld perturbations. The growth is slower than in 
standard cosmology and the growth exponent highly depends on a. For realistic values the growth 
rate changes from exp(4/3T) to between exp(1.28r) for a = 5 or exp(0.6r) for a = 3/2. This is a 
strong effect and it produces the unpleasant result that the structure formation should have taken 
a longer time for quintessence with exponential potential and the tracker solution. This seems to 
be a strong criteria to rule out or conflrm quintessential models. Further consequences should be 
investigated for this model and also for other potentials. 
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